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i Abstract 

The latest measurements of F2 at HERA allow for a combination of gluon and sea 
quark distributions at small x that is significantly different from those of existing parton 
sets. We perform a new global fit to deep-inelastic and related data. We find a gluon 
distribution which is larger for x < 0.01, and smaller for x ~ 0.1, and a flatter input sea 
quark distribution than those obtained in our most recent global analysis. The new fit also 
gives a s (Mz) = 0.114. We study other experimental information available for the gluon 
including, in particular, the constraints coming from fixed-target and collider prompt 7 
production data. 
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The advent of the high-energy electron-proton collider, HERA, has enabled the partonic 
structure of the proton to be investigated in the small x regime, x < 10~ 3 . The first measure- 
ments showed that the proton structure function F 2 (x, Q 2 ) rose dramatically with decreasing 
x. When these HERA data were included in global analyses of deep-inelastic and related 
data, the sea quark and gluon distributions were typically found to have the small x behaviour 
xS,xg ~ x~ ' 3 . To be precise, in the most recent MRS analysis the sea quark and gluon 
distributions are parametrized in the form 

xS = A s x- Xs (l + e s ^ + ls x)(l-x) r,s (1) 
xg = A g x' Xa {l + e 9V /x + 7 r r)(l - x) Va (2) 

with e g = and \$ = X g = A. At that time there were insufficient experimental constraints on 
the gluon to justify the introduction of the parameter e g in xg, or to determine the exponent 
X g independent of that of the sea-quark distribution S. The global fit gave A = 0.3 ± 0.1 
where the value was strongly correlated to eg. An independent global analysis by the CTEQ 
collaboration gave similar results with a common A of 0.286. Thus the small-x behaviour 
of the sea-quark and gluon distributions are closely linked - motivated either by BFKL gluon 
dynamics with the sea quarks driven by the g — > qq transition, or by GLAP evolution from a 
low starting Q 2 scale. 

The most recent measurements of F 2 (x,Q 2 ) presented by the ZEUS [§] and HI || collab- 
orations have more data points at low x and improved precision compared to those presented 
earlier. The new measurements are shown in Fig. 1, together with their description by MRS (A) 
partons (dash-dot curves). At very small x we see that the MRS (A) values lie above the 
new data. The same is true of the new CTEQ3 partons @ and, to a lesser extent, of the 
latest set of GRV(94) partons (dotted curves). We therefore repeat the global analysis of 
Ref. |l] with the new HERA data included. We allow the normalization of the ZEUS |4| and 
HI H data sets to vary within their quoted uncertainties (±3.5% and ±4.5% respectively). 
The resulting partons^, which we denote A', correspond to the following parametrization of the 
starting distributions at Qq = 4 GeV 2 



xu v = A t x a559 (l-0.54 v ^±4.65x)(l-£) 3 - 96 (3) 

xd v = A d x a335 (l±6.80 v / x± 1.93x)(l -x) 4 - 46 (4) 

xS = 0.956 x" 017 (l-2.55 v / x±11.2a;)(l-a;) 9 - 63 (5) 

xg = A g x-° A7 (l - l.90y/x + 4.07a;) (1 -a;) 5 ' 33 . (6) 



The sea S = 2(u+d + s + c) has a flavour structure similar to that of MRS (A), with parameters 
A A = 0.045, 7a = 49.6, r] A = 0.3, 5 = 0.02 and m 2 = 2.7 GeV 2 in the notation of Eqs. (7) 
and (8) of Ref. [J]. The QCD scale parameter is found to be A^(nf = 4) = 231 MeV, which 
corresponds to a s (M|) = 0.113. The parameters A u and Ad are determined by the flavour sum 
rules, giving 2.26 and 0.279 respectively, and A g = 1.94 is determined by the momentum sum 
rule. 

The A and A' partons differ only in the small-x region, since the data have changed only in 
this region. Indeed the main effect of the new HERA data is to change the small-x behaviour 

1 The FORTRAN code for the A' and G sets are available by electronic mail from W.J.Stirling@durham.ac.uk 
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A = \ s = X g = 0.3 of MRS (A) to A = 0.17 of MRS(A'). Besides the change in A, in the A' 
analysis we allowed e g to be a free parameter, which also helped to significantly improve the fit. 
The dashed curves of Fig. 1 show how the description of the small-x data is improved. However 
it is evident that the HERA data suggest stronger scaling violations than the new A' fit. Now, 
for x < 0.01, the slope 

mM2) - a - m £4 r d M-\ p„ (-) » 9 (»,q\ (7) 



<9iogQ 2 7T y qJx y \yJ \y 



is virtually a direct indicator of the magnitude of the gluon distribution. Indeed both the ZEUS 
and HI collaborations |7j have used their respective data to determine the slope, and hence 
the gluon distribution, and find values for the gluon in excess of MRS (A'), more consistent 
in fact with the steeply rising gluon distributions of, for example, MRS(D'_) ||. On the one 
hand the new A' fit to F 2 gives a flatter sea (and gluon) and yet on the other hand the slope 
dF 2 (x, Q 2 )/dlogQ 2 seems to prefer a steeper gluon. 

The situation is well illustrated by Fig. 2. The upper plot for dF 2 (x, Q 2 )/dlogQ 2 is an 
indicator of the gluon, see Eq. (0), whereas the lower plot is dominated by the sea quark 
distribution. The new HERA data show that the dashed A' curve, obtained assuming A g = As, 
is too flat (steep) in the upper (lower) plot. We conclude that the small-x data are now 
sufficiently precise that their description requires partons with A^ 7^ A5. We therefore repeat 
the global analysis with both A 9 and As free and find at Ql = 4 GeV 2 

xu v = A u x a593 (l-0.76 v / x + 4.20x)(l-x) 3 - 96 (8) 

xd v = A d x°- 335 (1 + 8.63 v / x + 0.32x)(l - x) AA1 (9) 

xS = 1.74x- a067 (l-3.45v^+10.3x)(l-x) iai (10) 

xg = A g aT°- 301 (l -4.14^x4- 10.1x)(l - x) 6 ' 06 , (11) 

with Ag = 1.51; A& = 0.043 and 7a = 64.9, the remaining parameters being unchanged from 
A'. The QCD scale parameter increases to 

A M g(n / = 4) = 255MeV, (12) 

which corresponds to a s (M|) = 0.114. The resulting set of partons 1 is denoted by MRS(G), and 
gives rise to the continuous curves in Figs. 1 and 2. We see that the steeper gluon (A g = 0.30) 
and the flatter sea (As = 0.07) significantly improve the description of the small-x HERA data. 
Due to the negative ^ g \fx term in xg, the effective value of X g is larger than the value shown. 
To be precise, if we approximate Eq. ([TTj) by xg = Ax~ x over the interval 10 -4 < x < 5 x 10~ 3 
then we find A = 0.355. Note that the data shown in Fig. 2(b) are the result of extrapolating 
measured data from higher values of Q 2 . Even if future low Q 2 measurements continue to show 
a rise for F 2 at small x the slopes extracted from Fig. 2(a) would still allow a steep gluon 
solution. 

The new HERA measurements of F 2 provide for the first time a reliable estimate of the 
gluon at small x. We can immediately see the trends. If it is indeed true that the HERA data 
prefer an enhancement of the A' gluon for x ~ 10 -3 [0 then there are implications for other 
gluon-dominated reactions, in particular for prompt 7 production. The fixed-target prompt 7 
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data |], |H pin the gluon at x ~ 0.4. This information, together with the momentum sum rule, 
implies (i) a decrease of the gluon at medium x ( ~ 0.1) and (ii) a steeper slope of the gluon 
from x ~ 10 _1 to x ~ 10 -3 . As a consequence we find (see below) that the gluon obtained 
from our new global analysis G does indeed change the description of the pp collider data on 
prompt photon production at small xt = 1p 1 T j\fs. 

Another important consequence of the reduction of the gluon in the interval x ~ 0.1 — 0.2 
concerns the prediction for the value of a s . At these x values the gluon still gives a significant 
(positive) contribution to <9F 2 /<91ogQ 2 . Since precise fixed-target measurements of F 2 exist in 
this x region (see for example Figs. 6 and 8 of Ref. Jl[]), the reduction in the gluon must be 
compensated by an increase in the value of a s . This brings the deep-inelastic determination of 
a s = 0.114 more in line with the world average value a s = 0.117 [{TCfl . 

Figure 3 compares the A, A' and G partons at Q 2 = 20 GeV 2 . We see that the quark 
distributions are essentially identical for x > 0.05. Indeed the A' and G description of the 
fixed-target deep-inelastic scattering data, the rapidity asymmetry and the asymmetry 
of Drell-Yan production in pp and pn collisions is as good as that obtained by MRS (A), see 
Ref. [I]]. However from Fig. 3 we see a sizeable difference in the A/ A' and G small-x gluon 
distributions. Fig. 4 highlights this difference, as well as indicating the x range of the various 
experimental constraints on the gluon. Comparing the G and A' gluons, we see (i) that the 
new HERA measurements of dF 2 /dlogQ 2 (Fig. 2) lead to an enhancement of the gluon for 
x < 0.01, (ii) that the fixed-target prompt photon data require the gluon to be unchanged 
for 0.35 < x < 0.55, see Fig. 5, and (iii) that as a consequence the G gluon is reduced in the 
intermediate interval 0.02 < x < 0.2. It is also interesting to compare the fraction of the gluon's 
momentum in various intervals of x for Q 2 = 20 GeV 2 for the two sets of partons: 



set 


< X < lO" 3 


10~ 3 < x < 10- 2 


10~ 2 < x < 10- 1 


lO^ 1 < x < 1 


MRS (A') 


3.5% 


13.4% 


46.1% 


37.0% 


MRS(G) 


5.9% 


16.0% 


42.3% 


35.7% 



Prompt photon production serves as a strong constraint on the gluon since it contributes at 
leading order. As mentioned above, the previous MRS analyses have used the WA70 fixed-target 
data H to pin down the gluon for x ~ 0.4. The procedure adopted was to perform a next- 
to-leading order fit with the renormalization and factorization scales optimized as described 
in Ref. \T1\. There now exists quite a range of prompt photon data which cover the interval 
0.01 < x < 0.6. Following Ref. fi~2H , we display these data in Fig. 5(a), together with their 
next-to-leading order description by the new G and A' parton sets with the renormalization 
and factorization scales set to fi — " " 
the R807 collaboration [IHf 



0.5pJ. In our plot we show the more recent ISR data from 
rather than the earlier R806 measurements. We note that these 
improved ISR data have a p^ T distribution whose shape is more in accord with QCD than those 
of R806. From Fig. 5 we see that, in general, the next-to-leading order QCD description of the 
prompt photon data is satisfactory, except that the p^ distributions observed by the CDF [13[ 
and UA2 |L5| collaborations are steeper than is predicted. Scale ambiguities could remedy a 
discrepancy in normalization but cannot significantly change the shape of a distribution. On 
the other hand the description of prompt photon production at high energies and small pj. is 
complicated by the presence of sizeable photon bremsstrahlung contributions. A detailed study 
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JTB[ of these photonic fragmentation effects (based on the full next-to-leading order calculations 
of Refs. [fT^I ) has been made, and as a result the QCD predictions are improved. However the 
description of the CDF and UA2 data shown in Fig. 5 already includes these next-to-leading 
order effects of photon fragment at ionP] and, although the shape is improved, it is still not 
satisfactory. It was hoped that the steeper MRS(G) gluon would have further improved the 
description of these data. We see from comparing the two plots of Fig. 5 that the improvement 



is, in fact, marginal. The most probable explanation [T2j of the residual discrepancy in the shape 
of the CDF and UA2 distributions at small xt is the presence of intrinsic kx arising from non- 
perturbative or 'soft' perturbative effects. Since no soft resummation calculations are available, 
we have modelled these effects by a two-dimensional transverse momentum convolution of the 
next-to-leading order QCD prediction for the jP T distribution with a Gaussian intrinsic k? form, 
exp(— /c^/Zcq), and varied the parameter fc 2 . to give the best fit to each set of prompt photon 
data. We find that an excellent description of the shape of the CDF and UA2 data sets can 
then be obtained if (k?) = y/nko/2 is taken to be 2.4 GeV in both cases. 

What other scattering processes could provide information on the gluon distribution at 
small x? In the production of bb pairs in high-energy hadron-hadron collisions the gluon enters 
(quadratically) at leading order via the subprocess gg — > bb. The prediction for the cross 
section is however very sensitive to the choice of ra& and the QCD scale, even at next-to-leading 
order, and reasonable variations of these parameters can change the normalization by a factor 



of 2 fL8| . Furthermore it is difficult experimentally to accurately reconstruct the 6-quark cross 
section from measurements of its decay products (£>, J/ip, leptons, ...). We have checked^ 
that for pp collisions at =1.8 TeV the differences in the A' and G gluons at small x only 
reveal themselves at pr{b) ^ 10 GeV/c where, unfortunately, there are essentially no CDF or 
DO data. At higher pr{b) the shape of the pr distribution measured by CDF |]18| is in reasonable 



agreement with the theoretical predictions, although the normalization is slightly larger than 
predicted with 'canonical' scale and mass values. 

Dijet production in pp collisions can also, in principle, probe the small x gluon p0| , pl| , [22| . 
For example, if two jets are produced with equal transverse momentum pt and pseudorapidity 
77 ^> 1 then x\ ~ 1, X2 ~ 2ptI\[s exp(— rf) < 1. A detailed calculation [[LlJ shows that at 
■s/s = 1.8 TeV the gluon can be probed in this way in the range 0.005 < x g < 0.05, as indicated 
in Fig. 4. However, at present the systematic errors are too large to allow any definite conclusion 
to be drawn. 

In a large region of phase space, inelastic J/ip production at HERA should be dominated 
by the 'colour singlet' process -yg — > J/ip g p3| . Here the gluon is probed at x ~ cMj/s 7P <C 1, 
where the constant c depends on the cuts used to define the inelastic events ||24j| . The theoretical 



calculation has recently been extended to next-to-leading order [25|. Forthcoming data from 
HERA, with for example 50 < y^s^ ^ 150 GeV and taking c = 3.4 f2H]| , will allow us to probe 
the interesting small x ~ 10~ 3 — 10~ 2 region, although it remains to be seen whether the 
precision will be high enough to discriminate between gluons like those of A' and G. 

Note that although there is a large overlap in the x region probed by pp dijet and HERA 
J/ip production, the Q 2 values are very different, being respectively Q 2 ~ 10 3 — 10 4 GeV 2 and 
Q 2 ~ 10 GeV 2 . Since the gluon distributions tend to become more similar as Q 2 increases, the 



2 We thank Werner Vogelsang for providing these predictions 
3 We thank Slawek Tkaczyk for performing these calculations 
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J/ip cross section at lower Q 2 provides greater discrimination. For example, the ratio of the 
G and A' gluons at x = 10" 3 is 1.71 at Q 2 = 5 GeV 2 and only 1.11 at Q 2 = 5 x 10 3 GeV 2 . 
The difference at the former Q 2 value is indicative of the difference in the inelastic J /if) cross 
section predictions at high y/s lp at HERA. 

In conclusion, we have seen that standard sets of next-to-leading order parton distributions 
like MRS (A), CTEQ3 and GRV(94) give a good overall description of the new HERA structure 
function data. However at the very smallest x values there is evidence that a new effect is 
becoming apparent, namely that the small- x behaviour of the gluon and sea-quark distributions 
are not linked at Q 2 = 5 GeV 2 . If, as in set MRS (A'), we force the distributions to be the 
same, then the slope dF 2 (x,Q 2 )/d\ogQ 2 at small x is underestimated. Our new MRS(G) set 
accommodates this by allowing the small-x quarks and gluons to have a different shape in x. 
In the notation of Eqs. (1) and (2) we find that X 9 = 0.30 and \$ = 0.07, though the precise 
values are correlated to the parameters e g and eg. The situation is summarized in Fig. 6, where 
we show the theoretical predictions and experimental measurements of the structure function 
F 2 at a 'typical' small- a; value, x = 0.0004. The A' and G curves are labelled with their gluon 
and singlet quark effective A values at this x and the starting value Q 2 = 4 GeV 2 . 

The resulting MRS(G) 'singular' gluon and 'flat' sea quark distributions do not have a 
ready explanation in terms of either perturbative or non-perturbative QCD. On the one hand 
GLAP evolution from a low scale, such as performed by GRV [||] (see also |2f|), develops both 
a steep gluon and a steep sea quark distribution at small x. The former is evident in the data, 
the latter is not. Again Fig. 6 summarizes the situation. The GRV(94) curve has a similar 
slope to MRS(G) (both gluons are steep) but overestimates the data (the GRV quarks are 
steeper than the G quarks). On the other hand it might be argued that the leading log(l/x) 
resummation, encapsulated in the BFKL equation, is more appropriate at small x. A singular 
(unintegrated) gluon is obtained, as required by the data, but again the steepness is fed, in 
this case via the fc^-factorization theorem, directly into the sea quark distribution and hence 
into F 2 . Of course the application of next-to- leading order, leading- twist QCD in the HERA 
small-x regime may be too naive. There are some indications that higher-order perturbative 
corrections may be non-negligible and that these may affect the structure function and its Q 2 



evolution in different ways |2^| . It is not impossible, for example, that the x — > behaviour 
of the quark and gluon distributions are more similar when higher-order corrections are fully 
taken into account. Unfortunately the theoretical technology necessary to investigate this with 
any precision is not yet available. If, on the other hand, the form of the starting distributions at 
Qq were primarily of non-perturbative origin, the dominance of the forward quark-proton and 
gluon-proton scattering amplitudes by pomeron exchange would again have suggested A# ~ X g . 
In any case, more precise structure function measurements at small x and moderate Q 2 will be 
very useful in shedding further light on this issue. 
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Figure Captions 



[1] The (a) ZEUS j|] and (b) HI || measurements of F2(x, Q 2 ) in the small x regime compared 
with the predictions of the MRS (A) partons of Ref. JI[ (dash-dot curve, shown only for 
x < 5.6 x 10~ 4 ), the MRS (A') and MRS(G) partons of this work (dashed and solid curves 
respectively), and the GRV(94) partons of Ref. || (dotted curve). 

[2] (a) The slopes, b(x), determined from linear fits, F 2 = a + blogQ 2 , to the HERA data 
and (b) the values of F 2 (x, Q 2 — 5 GeV 2 ) obtained by extrapolating the linear fits back to 
Q 2 = 5 GeV 2 . Also shown are the predictions of the MRS (A') (dashed curve) and MRS(G) 
(solid curve) partons of this work, together with the description given by the recent 
GRV(94) partons of Ref. (dotted curve). The slopes are calculated at Q 2 = 20 GeV 2 . 

[3] MRS(G) partons compared to MRS (A') and also to the earlier MRS (A) |jj partons, at 
Q 2 = 20 GeV 2 . For clarity, the gluons are shown for x < 0.005 only. 

[4] The MRS (A') and MRS(G) gluon distributions at Q 2 = 20 GeV 2 . Also shown are the x 
intervals in which the gluon is constrained by the various sets of data. 



[5] The description of WA70 @, UA6 [|, R807 @, UA2 |T§ and CDF @ prompt photon 
data by the MRS (A') (lower figure) and MRS(G) (upper figure) partons. The predictions 
for the pp collider cross sections are calculated using the next-to-leading order photonic 



fragmentation program of Ref. |L6|, and for the other cross sections using the next-to- 



leading order program of Aurenche et al. [T? 



[6] The description of the ZEUS [|J and HI || measurements of F 2 (x,Q 2 ) near x = 0.0004 
compared with the MRS (A') and MRS(G) partons of this work (dashed and solid curves 
respectively), and the GRV(94) partons of Ref. 0] (dotted curve). The effective A values 
for the gluon and singlet quark distributions at this value of x and Q 2 = 4 GeV 2 are 
indicated. 
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Zeus (xl.02) comparison at small x 
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